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Abstract 

The ideas of instantaneous amplitude and phase are 
well understood for signals with real- valued samples, 
based on the analytic signal which is a complex sig- 
nal with one-sided Fourier transform. We extend 
these ideas to signals with complex-valued samples, 
using a quaternion-valued equivalent of the analytic 
signal obtained from a one-sided quaternion Fourier 
transform which we refer to as the hypercomplex rep- 
resentation of the complex signal. We present the 
necessary properties of the quaternion Fourier trans- 
form, particularly its symmetries in the frequency do- 
main and formulae for convolution and the quater- 
nion Fourier transform of the Hilbert transform. The 
hypercomplex representation may be interpreted as 
an ordered pair of complex signals or as a quater- 
nion signal. We discuss its derivation and proper- 
ties and show that its quaternion Fourier transform 
is one-sided. It is shown how to derive from the hy- 
percomplex representation a complex envelope and a 
phase. 

A classical result in the case of real signals is that 
an amplitude modulated signal may be analysed into 
its envelope and carrier using the analytic signal pro- 
vided that the modulating signal has frequency con- 
tent not overlapping with that of the carrier. We 
show that this idea extends to the complex case, 
provided that the complex signal modulates an or- 
thonormal complex exponential. Orthonormal com- 
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plex modulation can be represented mathematically 
by a polar representation of quaternions previously 
derived by the authors. As in the classical case, there 
is a restriction of non-overlapping frequency content 
between the modulating complex signal and the or- 
thonormal complex exponential. We show that, un- 
der these conditions, modulation in the time domain 
is equivalent to a frequency shift in the quaternion 
Fourier domain. Examples are presented to demon- 
strate these concepts. 

1 Introduction 

The instantaneous amplitude and phase of a real- 
valued signal have been understood since 1948 from 
the work of Ville [l] and Gabor [2] based on the 
analytic signal. A critical discussion of instanta- 
neous amplitude and phase was given by Picinbono 
in 1997, particularly with reference to amplitude and 
phase modulation [3j §V]. The analytic signal and 
its associated modulation concepts can be simply de- 
scribed, even though a full theoretical treatment is 
quite deep. Given a real-valued signal f(t), the cor- 
responding analytic signal a(t) is a complex signal 
with real part equal to f(t) and imaginary part or- 
thogonal to f(t). The imaginary part is sometimes 
known as the quadrature signal - in the case where 
f(t) is a sinusoid, the imaginary part of the analytic 
signal is in quadrature, that is with a phase differ- 
ence of — tt/2. The orthogonal signal is related to f(t) 
by the Hilbert transform [4j[5]. The analytic signal 
has the interesting property that its modulus \a(t) \ is 
an envelope of the signal f(t). The envelope is also 
known as the instantaneous amplitude. Thus if f(t) is 
an amplitude-modulated sinusoid, the envelope |a(f)|, 
subject to certain conditions on the frequency con- 
tent, is the modulating signal. The argument of the 
analytic signal, Z a(t) is known as the instantaneous 
phase and its derivative is known as the instanta- 
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neous frequency. The analytic signal has a further 
very interesting property: it has a one-sided Fourier 
transform. Thus a simple method for constructing 
the analytic signal (algebraically or numerically) is 
to compute the Fourier transform of f(t), multiply 
the Fourier transform by a unit step which is zero for 
negative frequencies, and then construct the inverse 
Fourier transform of the resulting spectrum. 

In this paper we extend the above ideas to sig- 
nals with complex- valued samples. We show that 
if a quaternion Fourier transform is computed from 
a complex-valued signal, and negative frequencies in 
the frequency domain representation are suppressed, 
the inverse quaternion Fourier transform gives an 
equivalent of the analytic signal, which we refer to 
in this paper as a hypercomplex representation of the 
complex signal. Just as the classical analytic signal 
represents a real-valued signal by a complex signal 
(that is with a pair of real signals), the hypercomplex 
representation presented here represents a complex- 
valued signal by a quaternion signal (that is with a 
pair of complex signals, based on the Cay ley- Dickson 
form of a quaternion). 

We consider signals with complex-valued samples 
in a most general sense, without any restrictions 
on the relationship between the real and imaginary 
parts. However, we exclude the special case of an 
analytic signal (where the real and imaginary parts 
are orthogonal) because it does not possess an in- 
teresting hypercomplex representation: the two ad- 
ditional quaternion components of the hypercomplex 
representation are simply copies of the two original 
components. The case where the real and imaginary 
parts are correlated in some statistical sense is called 
improper, whereas an analytic signal is proper. It has 
been customary to handle such signals by using an 
augmented representation consisting of a vector con- 
taining the signal and its conjugate IB]. Recently, 
Lilly and Olhede [7] have published a paper on bi- 
variate analytic signal concepts. Their approach is 
linked to a specific signal model, the modulated ellip- 
tical signal, which they illustrate with the example of 
a drifting oceanographic float, and they utilise pairs 
of analytic signals derived from the pair of real signals 
representing the time-varying positional coordinates 
of a drifting float. In this paper a different approach, 
using hypercomplex representation, is considered be- 
cause it is more suited to the analysis of modulation, 
and also because we have at our disposal the existing 
tools of quaternion Fourier transforms, which provide 
a natural fit with modulation of a complex signal by 



a complex signal, for reasons which should be evi- 
dent later in the paper. Note that, in essence, the 
approach proposed in this paper is inspired by the 
work of Bulow 8,9 and Labunets flOl who consid- 



and Labunets 
ered Clifford Fourier transforms (quaternion in the 
case of Biilow) to process multidimensional signals. 
While Biilow considered 2D signals (images) , here we 
consider ID signals with 2D samples (complex val- 
ued samples). We are not focusing on any particular 
application in this paper, rather we propose a new 
approach to the processing of complex valued signals 
(including proper and improper signals, even though 
proper signals are well-described using the classical 
complex Fourier transform — they are a special case.) 
Schreier and Scharf's book gives examples of appli- 
cations [6]. 

The idea of amplitude modulation is extended to 
the case where the modulating signal is complex, and 
the 'carrier' is an orthonormal complex exponential, 
that is a complex exponential which is in a com- 
plex plane orthogonal to the complex plane of the 
modulating signal - this concept is easily realised 
using quaternions, since quaternion algebra has a 
4-dimensional basis with three mutually orthogonal 
imaginary units. Subject to the same restrictions on 
frequency content as in the real- valued case, ampli- 
tude modulation in the time domain corresponds to 
a frequency shift in the Fourier domain. Further, a 
complex envelope of the hypercomplex representation 
can be defined which recovers the complex modulat- 
ing signal. This is the complex instantaneous ampli- 
tude. Further, the phase of the orthonormal complex 
exponential may be recovered. This is the instanta- 
neous phase, and its derivative is the instantaneous 
frequency. In the case of modulation of a complex ex- 
ponential the instantaneous frequency is that of the 
complex exponential 'carrier'. Just as the classical 
complex analytic signal contains both the original 
real signal (in the real part) and a real orthogonal 
signal (in the imaginary part), our hypercomplex sig- 
nal representation contains two complex signals: the 
original signal and an orthogonal signal. 

We have previously published partial results on 
this topic 11 ■ 13 as we developed the ideas presented 



here. The idea of orthonormal complex modulation 
has not previously been presented. 

The sequence of topics in the rest of the paper is 
as follows. Section [2] presents the quaternion algebra 
and quaternion Fourier transform concepts used to 
construct the hypercomplex representation of a com- 
plex signal using a one-sided quaternion spectrum. 
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Section [3] discusses the hypercomplex representation 
of a complex signal. Section [4] discusses the instanta- 
neous amplitude and frequency concept for a complex 
signal through its hypercomplex representation. Sec- 
tion [5] presents some signal examples to illustrate the 
ideas presented in the paper. 

2 Quaternion Fourier 
Transform 

In this section, we present the definition and prop- 
erties of the Quaternion Fourier transform (QFT) of 
complex valued signals. Before introducing the main 
definitions, we give some prerequisites on quater- 
nions. Interested readers may refer to 14 15 for a 



more detailed presentation of quaternions. 



2.1 Preliminary remarks 

Quaternions were discovered by Sir W. R. Hamilton 



in 1843 14 . Quaternions are 4D hypercomplex num- 



bers that form a noncommutative division ring de- 
noted H. A quaternion q 6 H has a Cartesian form: 
q = qo + qii + q 2 j + q^k, with q , qi, q 2 , q 3 G K and 
i,j,k roots of —1 satisfying i 2 = j = k 2 = ijk = 
— 1. The scalar part of q is a: S(q) = qo. The vector 
part of q is: V(q) = q — S(q). We will also make use 
of the real and imaginary parts using the following 
notation: 



q = dt(q), qi=%{q), 



<72 = %(<?), 



93 = 3 fe (<?) 

(1) 

Quaternion multiplication is not commutative, so 
that in general qp 7^ pq for p, q G H. The conju- 
gate of q is q = qo — q\i — q 2 j — q 3 k. The norm of q is 

I Ml = M 2 = (lo +qf + Q2+ ?!) = QQ- A quaternion 
with qo = is called pure. \q\ is the modulus of q. If 
\q\ = 1, then q is called a unit quaternion. The in- 
verse of q is q^ 1 = q/\\q\\. Pure unit quaternions are 
special quaternions, among which are i, j and k. To- 
gether with the identity of H, they form a quaternion 
basis: {l,i,j,k}. 

In H, in addition to the conjugation, there exist 
anti-involutions which play an important role, as no- 
ticed by Biilow et al. 9 . In the special case where 
{l,i,j,k} is chosen as the basis for H these involu- 
tions are, for a given quaternion q, the following 



= -iqi, q d = -jqj, q k = -kqk. 



(2) 



Quaternions can also be viewed as complex numbers 
with complex components, i.e. one can write q = 



z \ + z 2j in the basis {l,i,j,k} with z%,z 2 € C l , i.e. 
z a = $t(z a ) + i^s(z a ) for a — 1,2. This is called the 
Cayley-Dickson form of q. 

Among the possible representations of q, we will 
make use in this paper of one that was recently in- 
troduced: the polar Cayley-Dickson form 16 . Any 



quaternion q also has a unique polar Cayley-Dickson 
form given by: 

q = A q exp(B q j) = (a + a x i) exp((6 + M)j) (3) 

where A q = aa + a\i is the complex modulus of q and 
B q = bo + bii its complex phase. As explained in [161, 
given a degenerate quaternion p = (c + di)j , then its 
exponential e p is given by: 



cos \p\ 



P_ 

\P\ 



sin |p| = a + f3j + 8k (4) 



leading to a = cos\p\, j3 — {c/\p\) sin \p\ and S — 
{d/\p\) sin with \p\ = \Jc 2 + d? . Recall also that, 
for any quaternion p, \e p \ — 1. Then, for any quater- 
nion q given in its polar form as in equation Q: 



cos \B„ 



■ 7 '|^| sin l jB '?l + fc j|rj sm l- B 9l ( 5 ) 



where \B q \ = \Jb\ + b\. As a consequence, using the 
right hand side expression in ([3]), given a quaternion 
q = qo+qii+q 2 j +q^k, then the complex components 
of its polar Cayley-Dickson form, i.e. A q and B q , are 
given as follows: 



B q = 



cos(^/ 'ql+qf) 



1 <il+1l 



(6) 



with K — ^/q 2 + q^ arctan \^\/q 2 + q 2 J ■ In Section 
we will make use of the complex modulus A q for tne 
interpretation of the hypercomplex representation. 

2.2 Quaternion Fourier transform 

In this paper, we use a ID version of the (right) QFT 
first defined in discrete-time form in [IT] . Thus, it 
is necessary to specify the axis (a pure unit quater- 
nion) of the transform. So, we will denote by QFT^ 
a QFT with transformation axis fi. In order to work 
with the classical quaternion basis {1, i, j, k} and for 
reasons that will be explained later, we will use j as 
the transform axis. The only restriction on the trans- 
form axis is that it must be orthogonal to the original 
basis of the signal (here {l,i}). We now present the 
definition and some properties of the transform used 
here. 
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Definition 1. Given a complex valued signal z(t) = 
z r (t)+iZi(t) that belongs to L 1 {R,C)nL 2 (R,C), then 
the quaternion valued function denoted Zj{v), and 
given as: 



+00 



Zj(u) = QFTj [z(t)] = J z(t)e-~^ vt At (7) 

— OO 

and with Zj{v) G L 1 (R,H) n L 2 (]R,II) is called the 
Quaternion Fourier Transform QFTj with respect to 
axis j of z(t). 

The inverse transform is denned as follows. 

Definition 2. Given a quaternion Fourier transform 
Zj{v) (with respect to axis j), then: 

+00 

z{t) = IQFTj [Zj(v)] = [ Zj{i>)e j27!Vt dv (8) 



is the Inverse Quaternion Fourier Transform IQFTj 
with respect to j of Zj(y). 

2.3 Properties of the QFT 

In order to provide insight into the advantages of us- 
ing the QFT to analyze complex signals, we present 
a non-exhaustive list of the QFT properties. 

Property 1. (Transform symmetry) Given a com- 
plex signal z(t) = z r (t) + izi(t) and its quaternion 
Fourier transform denoted by Zj{v), then the follow- 
ing properties hold: 

• The even part of z r (t) is in tft(Zj{i>)) 

• The odd part of z r (t) is in ^Sj(Zj(v)) 

• The even part of Zi{t) is in Qi(Zj(y)) 

• The odd part of Zi(t) is in ^Sk(Zj(y)) 

Proof. Expand (|7|) into real and imaginary parts with 
respect to i, and expand the quaternion exponential 
into cosine and sine components: 



%{Zj{y)) 







/ 














/ 










%(ZjM) 



% k {Zj(v)) 



Figure 1: Symmetry of the four components of the 
quaternion valued Fourier transform Zj(v) of a com- 
plex signal z(t). 



from which the stated properties are evident. 



□ 



It is a known fact that the classical Fourier trans- 
form of a real signal possesses Hermitian symmetry. 
It is also known that the classical Fourier transform 
of a complex signal possesses no symmetry at all. The 
QFT of a complex signal has symmetry according to 
the following properties: 

Property 2. (i-involution reversal) Given a complex 
valued signal z(t) and its quaternion Fourier trans- 
form denoted QFTj[z(t)] — Zj(y), then the following 
property holds: 



+00 



Zj(-v) = -iZj(v)i = Zj(v) 



(9) 



[z r (t) + izi(t)] [cos(27rz4) — j sin(27wt)] dt 



+00 +00 

z r (t) cos(2-Kist)dt — j J z r (t) sm(2iri>t)dt 

—00 

+ 00 +00 
i / Ziit) cos(2irvt)dt — k / Zi(t) sin(27rz4)d£ 



Proof. This can be directly checked from Property [T] 

□ 

This symmetry property of the QFT of a complex 
signal is illustrated in Figure [l] and is central to the 
justification of using the QFT to analyze a complex- 
valued signal carrying complementary but different 
information in its real and imaginary parts. Using 
the QFT, it is possible to keep the odd and even 
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parts of the real and imaginary parts of the signal in 
four different components in the transform domain. 
This idea was also the initial motivation of Biilow, 
Sommer and Felsberg when they developed the mono- 
genic signal for images [8,9 18 . Note that the use of 



hypercomplex Fourier transforms was originally in- 
troduced in 2D Nuclear Magnetic Resonance image 



analysis 19 20 



For use later 
property. 



in the paper, we mention another 



Property 3. (j -involution conjugation) Given a 
complex signal z(t) and its QFTj denoted Z j(v) , then 
the QFTj of its conjugate z* (t) is given by: 



QFT j [z*(t)} = -jZ j (v)j = Z j (v) J 



(10) 



Proof. This can be directly checked by calculation. 

□ 

An important issue for the forthcoming study of 
complex signals is the way modulation affects the 
QFT. Of special importance is ortho-complex mod- 
ulation. 

Property 4. (Ortho-complex modulation) Given a 
complex signal z(t) and its QFTj denoted Zj(v) 7 the 
QFTj of the ortho-complex modulation of z(t) by the 
quaternion-valued exponential exp(ji27Wgi) is 

QFTj[z{t)e j27UJat ] 



Zj(v-u ) 



(11) 



Proof. As e 



,j2irv t 



is an exponential with the same axis 



j as the QFTj , then direct calculation similar to the 
complex case can be performed as e^ a e^ b = e^ a+b \ 
The desired result is then directly obtained. □ 

Note that a similar frequency shift occurs simul- 
taneously in the four components of the quaternion 
valued Fourier transform. Note also that the axis of 
the modulation plays an important role as a modula- 
tion using axis i or k would not lead to a simple shift 
of the spectrum^] Also, it is obvious from this prop- 
erty and the convolution property (see equation |i"2|) 
fulfilled by the QFTj that QFTj[e^ Vot ] = 5{v-v Q ). 

Property 5. (Positive spectra modulation) Given 
a complex signal z{t) with a band-limited QFTj: 
Xj{v) = for \v\ > v m for a given v m . Then, for 
some vq > v m , the ortho-complex modulated signal 
z(t)e^ Vat has a right-sided QFTj. 



lr This is due to the fact that for two real numbers a and 
b, then e i a e kb ^ e ^ a+kb . This is in turn due to the fact that 
j and k do not commute. This is a speci al c ase of the well 
known Baker-Campbell-Hausdorff formula [2X| . 



-1*)) 



%{Zj(v - Vq)) 



^ k {Zj(v - v Q )) 



Figure 2: Four components of the quaternion val- 
ued Fourier transform Zj(y — v ) of the orthocom- 
plex modulated signal z{t)e 32vVot and where the com- 
plex signal z(t) has a QFTj verifying \Zj(v) \ = for 



> vq. This spectrum is an example of ( 11 ). 



Proof. This results follows directly from property [4] 
and the band-limited property of z(t). □ 

This property is illustrated in Figure [2] 
Now, we introduce the convolution property in the 
case where one of the signals is real valued. Note 
that we do not give the general case as only the above 
mentioned one will be of use in the sequel. 

Property 6. ( Convolution) Given g and f such that: 
g : R -> C and f : R -> M. the QFTj of their convo- 
lution product [g * f)(t) is: 

QFTj Kg * /)(*)] = QFTj [g(t)} QFTj [f(t)] (12) 

Proof. Taking the QFT • of their convolution product 



one gets the result shown in ( 13 ) which is true thanks 



to the fact that f(t) commutes with the j exponential 
as f(t) is real valued. □ 

Note that the order matters in the product of the 



QFTjS in (12 1 as they are quaternion and complex 



valued (which do not commute in the general case). 
This expression is thus valid only when / is convolved 
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+00 +oo 

QFT j [(g*f)(t)}= j j g(T)f(t-T)dTe- 2 1™' t dt 



-oo — OO 
+00 



+00 



J g{T)e' 2j ^ T dT J f(t')e- 2:invt 'dt' 

— 00 —00 

QFTjig^QFTjifit)} 



(13) 



on the right. Convolution on the left would result in 
the appearance of a conjugation in the Fourier do- 
main. 

2.4 Hilbert transform 

The quaternion Fourier transform of the Hilbert 
transformer x(t) = 1/irt is — jsign(^), where j is 
the axis of the transform. This arises simply from 
the fact that using the exponential with axis j in the 
transform will lead to a real and 3j part for the QFTj 
only, which is a degenerate quaternion isomorphic to 
a complex number. Substituting x(t) into ([7]), we get: 



QFTj 



+OO 



" 1 " 




r e —j2wvt 




= 1) 


t 



dt 



and this is clearly isomorphic to the classical complex 
case. The solution in the classical case is — isign(^), 
and hence in the quaternion case must be as stated 
above. 

It is straightforward to see that, given an arbitrary 
real signal y(t), subject only to the constraint that 
its classical Hilbert transform H [y(t)] exists, then one 
can easily show that the classical Hilbert transform of 
the signal may be obtained using a quaternion Fourier 
transform as follows: 

H [y(t)} = IQFTj [-j sign^)^ {v)\ (14) 

where Yj(v) = QFTj [y(t)}. This result follows from 
the isomorphism between the quaternion and com- 
plex Fourier transforms when operating on a real sig- 
nal, and it may be seen to be the result of a con- 
volution between the signal y(t) and the quaternion 
Fourier transform of x(t) = 1/irt. Note that j and 
Yj{v) commute as a consequence of y(t) being real. 
We can now define the Hilbert transform of a complex 
signal. 

Definition 3. Consider a complex signal z(t) — 
z r (i) + izi{t) and its quaternion Fourier transform 



Zj (y) as defined in Definition [7J Then, the hyper- 
complex analogue of the Hilbert transform of z(t), is 
as follows: 

Hi W)] = Fj 1 [-j signal (u) (15) 

where the Hilbert transform is thought of as: 
Hj [z(t)] — p.v. (z*"\)> where the principal value 
(p.v.) is understood in its classical way. This result 
follows from ( |14[ ) and the linearity of the quaternion 
Fourier transform. To extract the imaginary part, the 
vector part of the quaternion signal must be multiplied 
by —i. An alternative is to take the scalar or inner 
product of the vector part with i. Note that j and 
Zj(y) anticommute because j is orthogonal to i, the 
axis of Zj{v). Therefore the ordering is not arbitrary, 
but changing it simply conjugates the result. 

Note that the signal z(t) is considered to be non- 
analytic in the classical (complex) sense, that is its 
real and imaginary parts are not orthogonal. How- 
ever, this definition is valid if z(t) is analytic, as it 
can be considered as a degenerate case of the more 
general non-analytic case. We end this section with 
a useful property of the right-sided QFTj . 

Property 7. (Positive spectrum quaternion signal) 
Consider a QFTj denoted by Zj(y). Then, Zj(y) is 
said to be right-sided if Zj(v) — 0, W < 0. Now, 
if a QFTj denoted Zj (y) is right-sided, then its In- 
verse Quaternion Fourier transform IQFTj [Zj(y)\ 
is a quaternion valued signal. 

Proof. The QFTj of a complex signal has the sym- 
metry described in Property [2] and illustrated in Fig- 
ure [T] Assuming these symmetries, a right-sided 
QFTj can only be obtained from a QFTj, denoted 
Zj(y), by the following linear combination: (1 + 
sign{v))Zj{i>) which simply cancels out the negative 
frequencies. Now, the IQFTj of sign(^) is, as previ- 
ously seen, equal to j/itt. The convolution of \jirt 
with z(t) leads to a complex signal with 5i and 3j 
which, when multiplied by j becomes a degenerate 
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quaternion signal with only the Qj- and 3fc parts non- 
zero. Now, the total IQFTj of (1 + sign(v)) Z j (v) 
consists of z(t) + j * z(t)j (t) which is thus made 
of a real part and three imaginary parts of axis i, j 
and k. □ 

In the next section, we make use of the QFT and 
Hilbert transform to build up a hypercomplex repre- 
sentation for complex valued signals. 



3 The hypercomplex signal 
representation 

We define a hypercomplex representation, denoted 
z(t), of the complex signal z(t) by a similar approach 
to that originally developed by Ville [j]. The z(t) 
signal is quaternion-valued. The following definitions 
give the details of the construction of this signal. 

Property 8. (Hypercomplex extension) Given a 
complex signal z(t), one can associate to it a unique 
canonical pair corresponding to a (complex) modulus 
and phase. The modulus and phase are uniquely de- 
fined through the hypercomplex representation of the 
complex signal, which is quaternion valued. 

Proof. Cancelling the negative frequencies of the 
QFT leads to a quaternion signal in the time domain. 
Then, any quaternion signal has a modulus and phase 
defined using its CD polar form. □ 

Definition 4. Given a complex valued signal z(t) 
that can be expressed as z(t) = z r (t) + izi(t), then 
the hypercomplex representation of z(t), denoted z(t) 
is given by: 



z(t) = z(t)+jH j [z(t)} 



(16) 



where T-Lj [z(t)] is the hypercomplex analogue of the 
Hilbert transform of z(t) defined in the preceding def- 
inition. The QFTj of this hypercomplex representa- 
tion of the complex signal, denoted Zj(v), is thus: 

z j{ v ) = Zj{v) - j 2 sign(i/)Z,-(z/) 
= [1 + sign(^)] Z» 
= 2U(v)Zj(v) 

where U (y) is the classical unit step function. 

This result is unique to the quaternion Fourier 
transform representation of the hypercomplex repre- 
sentation, which has a one-sided quaternion Fourier 



spectrum. This means that the hypercomplex repre- 
sentation may be constructed from a complex signal 
z(t) in exactly the same way that an analytic signal 
may be constructed from a real signal, by suppression 
of negative frequencies in the Fourier domain. The 
only difference is that in the hypercomplex represen- 
tation, a quaternion rather than a complex Fourier 
transform must be used, and of course the complex 
signal must be put in the form z(t) — z r (t) + izi(t) 
which, although a quaternion signal, is isomorphic to 
the original complex signal. 

Property 9. (Hypercomplex phasor form) Given a 
complex signal z(t), one can associate to it a unique 
canonical pair [p(t), 4>{t)\ consisting of a modulus and 
a phase. The modulus and phase are uniquely defined 
through the hypercomplex representation of the com- 
plex signal. 

Proof. Cancelling the negative frequencies of the 
QFTj of a complex signal leads to a quaternion sig- 
nal in the time domain, with similar frequency con- 
tent (the original QFTj of the complex signal can 
be recovered thanks to the symmetry properties of 
the quaternion Fourier transform). As any quater- 
nion signal has a unique modulus and phase defined 
through its CD polar form, those are uniquely asso- 
ciated to the original complex signal. □ 

A second important property of the hypercomplex 
representation of a complex signal z(t) is that it main- 
tains a separation between the different even and odd 
parts of the original signal. 

Property 10. (Orthogonal split) The original signal 
z{t) is the simplex part ^ Theorem 1], J23] § 13.1.3] 
of its corresponding hypercomplex representation z(t). 
The perplex part is the orthogonal or 'quadrature ' 
component, o{t). They are obtained by: 



o{t) = \(z+{t)~W)) 



(17) 
(18) 



Proof. This follows from (16). Writing this in full by 



substituting the orthogonal signal for Hj [z(t)]: 
z{t) = z(t) + jo(t) = z r (t) + izi(t) + jo r {t) - koi{t) 



and substituting this into equation (17), we get 



z(t) = 



z r (t) 

-i [Zr(t) - 



- izi(t) 
izi(t) ■ 



-jo r (t) - koi(t) 
jo r (t) - koi(t)} i 
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and since i and j are orthogonal unit pure quater- 
nions, ij = —ji: 

_ 1 / z r (t) + izi(t)+jo r (t) - k 0l (t)\ 
~ 2 \+ z r (t) + iziit) - jo r (t) + koi(t) J 

from which the first part of the result follows. Equa- 
tion |l8| ) differs only in the sign of the second term, 
and it is straightforward to see that if z(t) is substi- 
tuted, z(t) cancels out, leaving o(t). □ 

4 Instantaneous amplitude, 
phase and frequency 

Using the hypercomplex representation z(t) of a com- 
plex signal z(t), we now investigate the concept of 
instantaneous amplitude and phase for z(t). 

Recall that z{t) is quaternion valued, and thus it 
can be expressed using its polar Cayley-Dickson form 
given in ^ as: 



s(t) = p{ty 



(t)i 



(19) 



where p(t) is called the instantaneous amplitude of 
z(t) and <f>(t) is the instantaneous phase of z(t). The 



pair 



p(t),^(t) 



is canonical and associated to z(t). 

The one-to-one correspondence between z(t) and the 
canonical pair can be demonstrated using the same 
arguments used by Picinbono in 3 . Indeed, as pre- 
sented in the previous sections, z(t) has a one-sided 
spectrum. This ensures the spectral characterization 
(see [3j II. B]) of the canonical pair. In addition, the 
one-to-one correspondence is ensured because e J< ^*) 
is a phase signal with respect to the QFTj. Now, just 

as explained in [3J III], e 3 '^*' is a unimodular signal, 
and in order for it to be one-sided, some restriction;]^] 
apply to 4>(t). In short, the phase 4>{t) is of the form: 



4>{t) = 8 + 2iris t + (j) b {t) 



(20) 



where 6 is a constant, is is a given frequency and 
<pb(t) is a band- limited contribution with maximum 
frequency B > vq . In Section [5j we present examples 
of the former two terms. 

With all these considerations, one can state the 
following. 



2 Note that the considerations given in [3] for instantaneous 
phase, amplitude and frequency and their link with the an- 
alytic (one-sided) signal are based on the Bedrosian theorem 
which extends easily to the QFTj case due to the properties 
presented in Section [2] and [3] 



Definition 5. Given a complex signal z(t), its in- 
stantaneous amplitude p(t) and instantaneous phase 
<p(t) are obtained through the polar Cayley-Dickson 
form of its hypercomplex representation z(t). 

As noted above, through this definition, p(t), 4>(t) 

is a canonical pair associated to z(t). Note that, un- 
like in the classical case of the analytic signal asso- 
ciated to a real signal, the instantaneous amplitude 
of z(t) is complex valued. However, it represents the 
low-frequency part (non modulated) of z(t). The in- 
stantaneous frequency can now be defined similar to 
the classical case. 

Definition 6. The instantaneous frequency of the 
complex signal z(t) is defined as: 



fit) 



1 d(j)(t) 
2tt dt 



(21) 



where <j)(t) is the instantaneous phase of z(t). 



Like the classical case, this definition follows natu- 
rally from the Taylor series for </>(£) about to, namely, 



0(f) = 0(to) + (t - t ) 



d(f>(t ) 
dt 



R 



dt 



dt 



where R is small when t is close to to and we see 
that d ^l"^ has the role of frequency. Unlike the clas- 
sical case, the complex nature of z(t) requires that 
the oscillation direction evolves over time (as can be 
seen for example in Figure [3] and [9| . Information on 
the orientation of this oscillation, i.e. the osculating 
plane orientation, can be recovered from the instan- 
taneous amplitude p{t). More precisely, the normal 
to the osculating plane, denoted n(t), is: 



hit) = p(t) x p'(t) 



(22) 



where p'(t) — dp(t)/dt and where x is the cross prod- 
uct of vectors corresponding to the time- varying com- 
plex signals: [z r (t) , Zi(t) , t] is the vector associated 
with the complex signal z(t) — z r (t) + iz^t). From 
a geometric point of view, p'(t) can be considered as 
the velocity of z(t). 

Such a geometrical description of a complex sig- 
nal can be achieved with the use of concepts defined 
through the QFTj. Such an analysis could be devel- 
oped to higher dimensional signals using geometric 
algebra Fourier transforms. 
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5 Examples 

We present three examples in which a band-limited 
complex signal A(t) modulates an orthogonal com- 
plex exponential exp(£?(i)j). This can be represented 
using the polar Cayley-Dickson form of ^ as: 



q{t) = A(t)exp(B(t)j) 



(23) 



The result is a quaternion- valued signal. We can ex- 
tract from this, two complex signals using the Carte- 
sian Cayley-Dickson form: q(t) = z(t) + o(t)j, where 
z{t) = z r (t) + Zi(t)i and similarly for o(t). In what 
follows we take the first of these to be our modulated 
complex signal (the second is a type of 'quadrature' 
signal) . 

In the first example, the orthonormal complex ex- 
ponential has constant frequency, so that simply: 



Bi(t) = 2-Kv t 



(24) 



where Vq is an arbitrary constant value. This example 
can be generated in the frequency domain, since the 
modulation in the time domain corresponds to a shift 
in the frequency domain by vq. The complex signal 
is zi(t). 

In the second example, i?2(i) is a step function of 
the form: 



B 2 (t) = 



if 
if 
if 



t E [0,0.25] 
t E [0.25,0.75] 
t e [0.75, 1] 



(25) 



where vq and v\ are arbitrary constant values. This 
example is generated in the same manner as the first 
one. The complex signal is thus z%{t) for this exam- 
ple. 

In the third example, the frequency of the complex 
exponential sweeps linearly from an initial value to 
double the initial value and back to the initial value. 
That is: 



B 3 (t) = a(n T/2 *n T/2 )(;-T/2) 



(26) 



where II is the boxcar function with T — 0.5 and a 
a constant in the present case. This example is also 
generated in the time domain. The complex signal is 
denoted z 3 (t). 

We present in the sequel how the modulating com- 
plex signal, and the frequency of the carrier can be 
recovered from the quaternion hypercomplex repre- 
sentation of the complex signal using the instanta- 
neous amplitude (the 'envelope') and the instanta- 
neous phase, both of which are computed from the 




Figure 3: Band-limited signal z\{t) (blue) and enve- 
lope pi (t) (black) of the hypercomplex representation 
zi(t) associated to z\(t). 



Cayley-Dickson polar form of the quaternion-valued 
hypercomplex representation of the complex signals 
zi(t), 22 (i) and z 3 (t). 

Figure [3] shows a modulated complex signal z\ (t) 
(in blue) constructed by modulating an orthonormal 
complex exponential with a band-limited complex 
signal. The latter was constructed from a pseudo- 
random signal by filtering in the frequency domain to 
suppress all frequencies higher than an upper limit of 
16 cycles over the length of the signal. The orthocom- 
plex exponential in this case has constant frequency 
four times that of the highest frequency of the base- 
band complex signal, and the modulated signal was 
created by frequency shifting the quaternion Fourier 
transform of the baseband signal, as in (111. The 



complex envelope pi(t) of the hypercomplex repre- 
sentation of zi(t), i.e. z\(t), is displayed in black in 
Figure [3] 

Figure [4] shows the instantaneous phase <f>i(t) of the 
hypercomplex representation zi(t) of the signal Z\{t) 
displayed in Figure [3j This phase is obtained using 
the polar Cayley-Dickson phase of the hypercomplex 
representation z\(t). 

In the second example, the orthocomplex exponen- 
tial is again constant, but this time 'by parts', with 
two different constant values taken over the time pe- 
riod considered, as expressed in (25 1. The initial fre- 



quency of the exponential is 12.5 times that of the 
highest frequency of the baseband complex signal, 
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Figure 4: Phase 4>\(t) recovered from the polar 
Cayley-Dickson form of the hypercomplex represen- 
tation z\ (t) of the signal Z\ (t) displayed in Figure |3j 
showing the linearly increasing phase B\ (t) . 



Figure 6: Phase ^(i) recovered from the polar 
Cayley-Dickson form of the hypercomplex represen- 
tation 22 (t) of the signal Z2(t) in Figure [5] after un- 
wrapping. 




Figure 5: Band-limited signal 22(f) (blue) enve- 
lope f>2 (t) (black) of the hypercomplex representation 
22(f) associated to 22 (t). 



and it then doubles at time 0.25 before returning to 
its original value at time 0.75. Figure [5] shows the 
modulated signal in this case. The recovered enve- 
lope P2(t) is exactly as in the previous example, the 
frequency changes in the complex exponential makes 
no difference to the recovery of the envelope. The po- 
lar Cayley-Dickson phase 02 (f) is plotted in Figure^ 
after unwrapping. 

For this second example, we also present two 
zoomed views, Figures [7] and [8j taken for two different 
time intervals. On these two figures, in addition to 



the original signal z 2 (t) and the envelope /S 2 (f), some 
local geometric information is displayed. At several 
times, the vectors p' 2 (t) (blue) and n 2 (f) (red) are dis- 
played. One can identify that p' 2 (t) is the vector tan- 
gent to the envelope p%{t) of the signal z 2 (t), while 
the vector n 2 (£) is normal to the osculating plane. 
This local information gives access to the plane in 
which the signal is 'oscillating' over time. Note that 
this geometric information has been obtained from 
attributes derived from the hypercomplex represen- 
tation of the original signal 22 (f), i.e. through sim- 
ple spectral considerations. Such geometric quanti- 
ties are classically obtained from 3D curves using the 
Frenet-Serret formulae [24j Chap. I] for example. 

In the third example, the orthocomplex exponen- 
tial has a much higher frequency than the baseband 
signal (which is the same signal as in the first ex- 
ample). The initial frequency of the exponential is 
25 times that of the highest frequency of the base- 
band complex signal, and it sweeps linearly to dou- 
ble the initial frequency, and then back down to the 
initial value. The original signal 23 (t) together with 
the complex envelope p?, (f ) of the hypercomplex rep- 
resentation 23 (£) are displayed in Figure [9] In Fig- 



ure 
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we present the instantaneous frequency <fi' 3 (t) 
obtained from 23 (£). The linearity of the chirp fre- 
quency behaviour is recovered (singularity at time 0.5 
is due to the behaviour of the phase at this point) as 
expected. 
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Figure 7: Original signal 22 (i) (blue) and complex en- 
velope P2(t) (black) of the orthocomplex modulation 
[5] for the time period: t = 0.3 to t = 0.6. For dif- 
ferent times, the tangent vectors (blue), normal (red) 
and binormal (green) are displayed. The blue vec- 
tors represent p^t) for some given t, the tangent to 
the complex envelope. Green vectors represent fi2{i) 
(given in (22)), i.e. the normal to the instantaneous 



osculating plane of the signal 22 (i). 




Figure 8: Original signal Z2(t) (blue) and complex en- 
velope P2(t) (black) of the orthocomplex modulation 
[5] for the time period: t = 0.55 to t = 0.8. For dif- 
ferent times, the tangent vectors (blue), normal (red) 
and binormal (green) are displayed. The blue vec- 
tors represent p' 2 (t) for some given t, the tangent to 
the complex envelope. Green vectors represent n 2 (t) 
(given in (22)), i.e. the normal to the instantaneous 



Figure 9: Band-limited signal z%{t) (blue) enve- 
lope P3 (t) (black) of the hypercomplex representation 
zz{t) associated to 23(f). 




osculating plane of the signal z^{t). 



Figure 10: Instantaneous frequency 4>' 3 (t) obtained 
from the hypercomplex representation Zs(t) of the 
complex signal z 3 (t). 



Note that the notion of instantaneous frequency 
that we have introduced here for complex signals is 
indeed very similar to the notion of angular velocity 
that is well known in physics. The analogy is pos- 
sible thanks to the fact that the knowledge of n(t) 
gives access to the time- varying orientation of the 2D 
plane in which the signal is oscillating. 

Finally, a short-time quaternion Fourier transform 
(STQFT) of the modulated signal 23 (t) is shown in 
Figure |11| showing the frequency sweep of the ortho- 
complex exponential. The modulus of the STQFT 
is real valued and it shows the frequency content 
evolution over time for z$(t). The values for which 
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50 100 150 200 250 300 350 400 450 

Figure 11: Time-frequency depiction of the modu- 
lated signal z 3 (t) in Figure [9] The figure shows the 
modulus of a short-time quaternion Fourier transform 
with window length 128. 

the STQFT modulus vanishes represent the moment 
when the complex envelope vanishes. 

Notice that we have not illustrated the case where 
B q (t) is complex. This is because it is the modulus 
of B q (t) which is the phase in the orthocomplex ex- 
ponential, and a complex B q (t) would need to have a 
modulus representing phase. It is not clear what the 
argument would represent. A constant frequency is 
represented by a linearly increasing real phase, that 
is B q (t) would be a real- valued ramrj^J 

As can be seen from if B q (t) is real, that is 
B q (t) = b Q (t), we have the following form for the 
modulated signal, where A q (t) = ao(t) + ai(t)i: 

q(t) = A q (t) [cos |6o(t)| + sin \B (t)\j] 

= A q {t) cos |6o(t)| + A q (t) sin \B (t)\j 

This quaternion signal contains A q (t) modulated by 
a cosine in the first Cayley-Dickson component (the 
first two components of the quaternion), and A q (t) 
modulated by a sine in the second Cayley-Dickson 
component (the third and fourth components of the 
quaternion). The orthonormal formulation using 
quaternions achieves the following: 

• The 'carrier' is modeled by a complex exponen- 
tial, rather than a trigonometric function, giving 
the usual advantages of the notation. 

3 In fact it could be an imaginary- valued ramp with minimal 
change to the result. 



• Because the quadrature component is orthogonal 
to the in-phase component, the two are separated 
in the signal and in the Fourier domain. 

These advantages disappear if we replace j with i 
and revert to complex algebra, because then the in- 
phase and quadrature components become mixed in 
a single complex signal, rather than being kept sep- 
arate in the components of a quaternion signal. In 
complex algebra we can represent a complex expo- 
nential modulated by a real signal, or a real 'carrier', 
represented by a trigonometric function, modulated 
by a complex signal. We cannot represent a complex 
signal modulating a complex exponential. 

6 Conclusions 

We have shown in this paper how the classical con- 
cepts of instantaneous amplitude and phase may 
be extended to the case of complex signals using a 
quaternion hypercomplex representation of a com- 
plex signal. We have also shown that the classi- 
cal case of amplitude modulation can be extended 
to the complex case provided that we maintain an 
orthogonal separation between the complex plane of 
the modulating signal, and the complex plane of the 
complex exponential 'carrier', and maintain a separa- 
tion in frequency between the modulating signal and 
the 'carrier' exactly as in the classical case. The or- 
thogonal separation required is easily realised using 
a quaternion representation, and we can then com- 
pute a hypercomplex representation of the complex 
signal using a suitable quaternion Fourier transform. 
This hypercomplex representation of the signal pre- 
serves symmetries that would be lost using a complex 
Fourier transform. It is analogous to the classical an- 
alytic signal. 

We have shown that the quaternion polar repre- 
sentation permits us to recover from the hypercom- 
plex representation, both the instantaneous ampli- 
tude and the instantaneous phase. From the instan- 
taneous phase, it is also possible to recover the in- 
stantaneous frequency as well as the osculating plane 
that gives geometric information about the direction 
in which the signal is oscillating. 

As quaternions are isomorphic to the geometric 
Clifford algebra C/(0, 2) of ]R 2 , it is evident from our 
results that there are possibilities to extend the work 
presented here to the case of vector- valued signals 
where the samples are TV-dimensional, by using Clif- 
ford Fourier transforms in 2 N dimensions. 
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